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Abstract. In this paper we provide theoretical results that relate steady 
states of continuous and discrete models arising from biology. 



1. Introduction 

Mathematical modeling has been used successfully to understand the behaviour 
of many biochemical networks, [2 [Sj [191 US HI HHl [3 [3]. Mathematical models can 
be used as a framework to identify the key parts of a biological system, to perform 
in silico experiments and to understand the dynamical behaviour. 

There are mainly two mathematical frameworks used in modeling biological sys- 
tems. The continuous framework, usually constructed by coupling differential equa- 
tions, [2j [5l [191 [111; and the discrete framework, usually constructed by coupling 
multistate functions (e.g. Boolean functions) , [Tl [T8l l9l[3]. 

For many biological systems that have been modeled using both frameworks, it 
has been shown that both models have similar dynamical properties, [H [1] [TOl 131 
I13[ 118] : this raises the question of when and why this happens. Furthermore, it 
has been hypothesized that the dynamics of biochemical networks are constrained 
by the topological structure of the wiring diagram, [16] . This in turn suggests that 
mathematical models describing the same biochemical phenomena should have sim- 
ilar behaviour, even if they come from different frameworks. By providing answers 
to these questions we can gain substantial biological insight and understand how 
biology works at the system level. 

This problem has been studied by several authors, [H [20l [TOl [14]. Roughly 
speaking, the idea is that if a continuous model is "sigmoidal enough" and has the 
same qualitative features of a discrete model, then both models should have similar 
dynamics. In this paper we provide simple proofs about the relationship between 
steady states of continuous and discrete models satisfying these properties. Our 
results generalize previous results in [20l [H [TOl [M]. Furthermore, our proofs also 
show why the relationship between continuous and discrete models is likely to occur 
even when continuous models are not very sigmoidal. 

2. Preliminaries 

For a biochemical network with N species, a mathematical model is typically 
of the form x{t) = (xi(i), . . . , XAr(t)); where Xi{t) denotes the concentration of 
species i at time t. Given an initial condition x(0) = xq, it is of interest to study 
the behaviour of x{t) as t increases. There is usually no closed form of Xi{t) as a 
function of t (and xq), but rather equations that describe how x^'s depend on each 
other. 
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A system of ordinary differential equations (ODE) that models a biochemical 
network is typically of the form a;J(t) ~ fi{x{t)) — diXi(t) , where di is the elimination 
rate constant, corresponding to natural decay, and fi is a function that determines 
how the rate of change of Xi at time t depends on the value of the other species. 
We can represent this ODE model as x' ~ f{x) - Dx, where / : ^ and D 
is a diagonal matrix. 

A multistatc network or multivalued network (MN) that models a biochemical 
network is typically of the form Xi{t + 1) = fi{x{t)) (notice that time is discrete); 
where .t,; can take values in the finite set [0,7^,;] = {0, 1, . . . , m,} and fi : S = 
[0,mi] X |0, 7712] X ... X [0,mAr]] — >■ |0, m^] determines how the future value of 
species i depends on the current value of the others. We can represent this MN 
model as x{t + 1) = f{x{t)), where / : 5 — >■ 5. A multistate network is also 
called a finite dynamical system. Examples of multistate networks include Boolean 
networks, logical models and polynomial dynamical systems |17| . 

For any modeling framework, the wiring diagram or dependency graph is defined 
as the graph with N vertices, {1, . . . , N}^ and an edge from j to i if fi depends on 
Xj. The edge is assigned a positive (negative) sign if fi is increasing (decreasing) 
with respect to Xj . 

From now on we will use the notation x and / to refer to continuous entities; x 
and / will denote discrete entities. Also, N will denote the dimension of functions. 
II II will denote the Euclidean norm for vectors and the Frobenius norm for matrices. 
We denote closed intervals by [a,b], open intervals by ]a, &[ and half-open intervals 
by [a, b[ and ]a, b]. 

Our goal is to relate the steady states of an ODE and a MN; that is, fixed 
points of / and /. We assume (rescaling if necessary) that the ODE has the form 
x' = D{f{x) — x) where is a positive diagonal matrix and / : [0, 1]^ [0, 1]^. 
A MN is of the form f : S ^ Jlj^i lO, mjj S. 

Example 2.1. Our toy example consists of a 2-dimensional ODE, x' = f{x)~x, and 
a 2-dimensional MN, x{t + 1) = f{x(t))] where / : [0, 1]^ [0, 1]^ and / : [0, 2f ^ 
[0, 2p are given by: 
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where n = (ni, . . . , ng) is a parameter. Also, / is given as a truth table: 



Xl X2 I fi{x) f2{x) 





10 2 

2 1 2 

10 2 

112 2 

12 2 2 

2 2 1 

2 12 2 

2 2 2 2 



STEADY STATES OF CONTINUOUS AND DISCRETE MODELS 



3 



Notice that for our toy example D is the identity matrix. We wiU refer to the 
discrete values 0, 1 and 2, as low, medium and high, respectively. 

Figure [T] shows plots for /i and /2. We can sec that / and / have the same 
qualitative features. For example, when the value of both inputs is "low" , the value 
of both outputs is low; when the value of both inputs is high, value of both outputs 
is high. For a better comparison, heat maps for all functions arc shown in Figure 
Eland Figure 13 
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Figure 1. Graphs of functions /i (left) and J2 (right) from Ex- 
ample EHJ All exponents n^'s have been set equal to 10 for an 
easier comparison with the MN. 




Figure 2. Heat maps for function /i (left) and /i (right) from 
Example 12.11 L, M and H denote low, medium and high, respec- 
tively. 



Notice how the heat maps have the same qualitative features. The heat maps 
also show how to divide the domain of / to capture its qualitative behaviour. For 
example, the values of xi can be separated in three regions: [0,.3[ (low), ].3,.6[ 
(medium) and ].6,1] (high). Similarly, the values of X2 can be separated in three 
regions: [0,.4[ (low), ].4,.7[ (medium) and ].7,1] (high). 

Definition 2.2. Consider k = (fci, ^2, . . . , fc„j) e (0,1)'" such that < fci < 
^2 < . . . < fem < 1. For each element of |0, m\ — {0, 1, ... , to} we define [0]^ = 
[0,fci[, [l\k fc2[, • ■ • , [m - l]fe =]fcm_i,fc„[, [m\k =]fc„,l]. The subscript will 

be omitted if k is understood from the context. The numbers kiS will be called 
thresholds. 

Now consider fc' G (0,1)'"' for i — l,2,...,N. For each x S nr=ilOi'^«l 
define [x] = [xij^i x [x2\k'^ x . . . x [xjv]fc«. The sets [x\s will be called regions. 
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Figure 3. Heat maps for function /2 (left) and /2 (right) from 
Example 12.11 

For instance, from Example 12.11 we got the thresholds k-^ = (.3, .6) for the first 
variable and P (.4, .7) for the second variable. Then, [00] = [0, .3[x[0, A[, [01] = 
[0, .3[x].4, .7[, [02] = [0, .3[x].7, 1] and so on. AU regions are shown in Figure HJ 
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Figure 4. Regions for the domain of / in Example 12.11 



3. Results 

In this section we will prove that under certain conditions, there is a one-to-one 
correspondence between the steady states of ODEs and MNs. 
We first need the following lemmas. 

Lemma 3.1. Let f : [0, 1]^ — > [0, 1]^ be a continuous function. If K is a convex 
compact set such that f{K) C K , then the ODE x' = D{f{x) — x) has a steady 
state in K . 

Proof. It is a direct application of Brouwer's fixed point theorem. □ 

Lemma 3.2. Let f : [0, 1]^ — > [0, 1]^ he a differentiahle function. If K is a convex 
subset o/[0, 1]^, and \\f'{x)\\ < e < I for all x G K, then the ODE x' = D{f{x)-x) 
has at most one steady state in K . 

Proof. If x,y £ K are steady states of x' = D{f{x) — x), then f{x) = x and 
f{y) = y. Also, — y\\ = \\f{x) — f{y)\\ < e\\x — y\\; it follows that x = y. □ 

Lemma 3.3. Let / : [0,1]^ -> [0,1]^ be a differentiable function. If ||/'|| < 
on K , then any steady state of x' = D{f{x) — x) in K is asymptotically 

stable. 
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Proof. Suppose that x E K is a. steady state of x' = D{f{x) — x) and denote 
B = Df'ix). Notice that the Jacobian matrix of D{f — x) \s B — D. Suppose A is 
an eigenvalue oi B ^ D. Then, by the Gershgorin circle theorem, there is k such 
that I A — {Bkk ~ -Dfcfc)! < X^j^/c \-^kj ~ Dkj\- Since I? is a diagonal matrix, we obtain 
|A - {Bkk - Dkk)\ < Ej^k \Bkj\. It follows that |A + Dkk\ < E \Bkj\ < VN\\B\\. 

Then, \X + Dkk\ < Vn\\B\\ = ^\\f{x)D\\ < VN\\f'{x)\\\\D\\ < min{Aj < Dkk; 
thus we obtain |A + Dkk\ < Dkk- It follows that A must have negative real part. 
Therefore, x is an asymptotically stable steady state. □ 

Using the lemmas above we can easily prove the following theorem that relates 
steady states of continuous and discrete networks. Condition (1) states that the 
continuous and discrete functions have the same qualitative behaviour (for n large 
enough). Condition (2) states that the continuous function is "sigmoidal enough". 
We use n to denote a vector parameter and limits refer to all entries of n going to 
oo; for instance, in Example 12.11 n ~ (ni, . . . ^hq) and limits refer to — ^ oo for 
i = l,2,...,6. The parameter n controls how sigmoidal the function is. 

Theorem 3.4. Let (/") be a family of continuous functions from [0, 1]^ to itself 
and let f : S ^ S he a MN. Consider the following conditions: 

(1) /" ""''^> Lx G [./(a;)] on compact subsets of [x], for all x E S. 

(2) /" ""''^> on compact subsets of [x], for all a; G 5. 

Now, for each x Cz S consider a convex compact set C \x\ such that Lx is 
an interior point of Kf^^x) (with the topology inherited from [0, 1]^ ) and let K = 
UxesK.- 

Lf (1) holds, then for n large enough there is a one-to-one correspondence between 
steady states in K of the ODE x' — D(f"'{x) — x) and f. Furthermore, there is a 
steady state in Kx if and only if x is a steady state of f . Also, if x^^ G Kx is the 
steady state of f" , we have x" — > Lx. Additionally, if (2) holds, such steady states 
are unique and asymptotically stable. 

Proof. Suppose (1) holds. Since /" '""'^> Lx on Kx and Lx is an interior point 
of LC f(x) for all x G 5, there exists rig such that for all n > uq and for all x € S, 
f^{Kx) C Kf(^x)- If a; is a steady state of /, that is f{x) = x, then f^{Kx) ^ Kx', 
by Lemma [3.11 there exists a steady state of the ODE in Kx- On the other hand, 
if X is not a steady state of /, then Kx n f"'{Kx) C Kx n Kf(^x) ^ [A H [f{x)\ = {} 
and there cannot be a steady state in Kx- 

Now, for X € S such that f{x) = x, denote with x" a steady state of /" in Kx- 

Then, since /" ""'•^> Lx on Kx, we have that = /"(x") — Lx. 

Now suppose that (2) also holds. Consider e < "^^^^'^^ < 1 Since /" '""'^> 

on Kx, there exists ni > uq such that for all n > tii, ||/" || < e on I\. Then, by 
Lemma l3.2l and[ 3.31 the steady state must be unique and asymptotically stable. □ 

Example 3.5. It is straightforward to check that the function in Example l2 . 1 I satisfics 
the following {K denotes a compact set): 
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0, 0) e [00] 


= [/(oo)], 


[x] 


= [00] 


0, .9) e [02] 


= [/(oi)], 


[x] 


= [01] 


6,.9)e[12] 


= [/(02)], 


[x] 


= [02] 


8, 0) G [20] 


= umi 


[x] 


= [10] 


8, .9) e [22] 


= [/(II)], 


[x] 


= [11] 


8, .9) e [22] 


= [/(12)], 


[x] 


= [12] 


8,.5)G[21] 


= [/(20)], 


[x] 


= [20] 


8, .9) e [22] 


= [/(21)], 


[x] 


= [21] 


8, .9) e [22] 


= [/(22)], 


[x] 


= [22] 



Also, it is easy to see that /' !• (as n — > oo) on compact subsets of any region. 

Therefore, Theorem 13.41 guarantees that for n large enough we have a one to one 
correspondence between steady states and they are asymptotically stable. Since 
00 and 22 are the steady states of /, the ODE has two stable steady states. For 
n = (2,2,2,2,2,2) we obtain the steady states (0,0) G [00] and (.73, .78) € [22], 
shown in Figure [S] The phase portrait of / is constructed by placing an arrow from 
X E S to y G S ii f{x) = y; steady states are denoted by dots. 




Figure 5. Phase portrait of / (left) and / (right). Steady states 
of / are denoted by dots. 

It is important to mention that Theorem 13.41 generalizes previous results. For 
example, by restricting the theorem to piecewise-linear differential equations we 
obtain Theorem 1 in |14j : by restricting the theorem to Boolean networks and Hill 
functions we obtain Theorem 2 in [50]. 

The following corollary states that the set where we might not have a one-to-one 
correspondence can be made as small as possible. 

Corollary 3.6. Suppose that condition (1) from Theorem \3.4\ is satisfied and con- 
sider e > 0. There exists K C [0, 1]^ with /i([0, 1]^\ A") < e (n =Lebesgue measure) 
such that for n large enough there is a one-to-one correspondence between steady 
states in K of the continuous and discrete network. 
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Another implication of Theorem 13.41 is that the steady states of an ODE x' = 
D{f{x) — x) are located either near Lx or near the thresholds. This approach 
was used in [TUj to estimate the stable steady states of an ODE model for Th-ccU 
differentiation using as a starting point the steady states of a discrete model. Our 
results support this heuristic approach. 

It is important to mention that although Thcorem l3. 41 indicates that n has to be 
large, which would be meaningless for real parameters in biological regulation, the 
conditions for Lemma [3.11 13.21 13.31 can be satisfied in practice for low values of n. 
For instance, it turns out that the conclusion of Theorem 13.41 holds for values of rii 
as low as 2 (see Figure [5]). 

In other words, a continuous function can be sigmoidal enough for biological 
meaningful parameters. This can explain why many continuous and discrete mod- 
els of biological systems have similar behaviour [21 [1] [191 131 [IHl [12 . This also 
supports the conjecture that the dynamics of biological systems strongly depend 
on the "logic" of the regulation and not on the exact kinetic parameters, [HI [3]; in 
particular, the dynamical behaviour of biological systems is very robust to changes 
in the parameters. 



4. Application 

In this section we show how our results can be used to gain understanding on 
how the dynamics of a continuous model depends on the wiring diagram. Other 
applications of results relating discrete and continuous models have been shown in 

[2Ql[lo]. 

It is a well known fact that the topology of a network can constrain its dynamics. 
For example, it has been shown that positive feedback loops are responsible for 
multistationarity [T5l [TTl [T2] . 

We now will use our results to apply a theorem about MN to ODEs. First we 
need the following terminology. Consider a (signed directed) graph, G = {Vc,Ec) 
with vertices Vg and edges Eg- A positive feedback vertex set (PFVS) is a set 
P C Vg such that it intersects all positive feedback loops. In [3] and [T^ the 
authors proved the following theorem for Boolean networks and MN, respectively. 

Theorem 4.1. Let f : S = U^^iP^'^jl S be a MN and let P be a PFVS. 

Then, the number of steady states is bounded by Yiiepi^i Notice that in the 

Boolean case we have the bound 21^'. 

Combining our results and Theorem 14. II we easily obtain the following result. 

Theorem 4.2. Under the assumptions of Theorem \3.4\ and for n large enough, the 
number of steady states in K of x' = D{f"(x) — x) is at most Y[i£pi''^i + l); where 
P is a PFVS of the wiring diagram and the set [0, 1]^ \ K can be made as small as 
required. 

Example 4.3. Consider the differential equation x' — D{f{x) — x) where D is a 
positive diagonal matrix and / is given by: 

f? = 



/l = 




.5"! 


.5' 








.5"5 


.5' 


r— — ttk- 


/7 = 






.5' 





J5 



/8 = 



h = 




.5"4 


.5' 




h = 






.5' 




h = 






.5' 


'13+X-13 
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Figure 6. Wiring diagram of / in Example 13.51 



The wiring diagram is shown in Figure |6l Similarly to Example 13.51 it is not 
difficult to check that / has the same qualitative properties as a Boolean function 
(the thresholds are .5 in this case). Also, it is easy to check that the set P = {5} is 
a PFVS. Therefore, by Theorem 14. 1[ for n large enough wc have at most 21^' = 2 
stable steady states in K. 



The problem of relating continuous and discrete models has been studied by 
several authors, [5J HOI IIHI [SI- Previous results focused on piecewise linear and 
Boolean functions. We have shown that for sigmoidal ODEs, a steady state in the 
discrete model gives rise to a steady state in the continuous model; furthermore, this 
steady state is unique and asymptotically stable. Our results generalize previous 



One application of our results is the ability to extend the applicability of tools 
about discrete models to continuous models as shown in Section |4l The problem 
of relating network topology to dynamics has been studied extensively for discrete 
models; those results can give new insight on how network topology constrains the 
dynamics in continuous models. 

A natural question arising from our work is whether or not one can obtain similar 
results about periodic solutions. For some biological systems and in very partic- 
ular cases it has been shown that continuous and discrete models produce similar 
periodic behaviour [SI [71 [H [TJ [Tl] . General results in this direction would increase 
our understanding of the relationship of continuous and discrete models, how the 
dynamical properties are constrained by topological features of the wiring diagram 
and how biology works at the system level. This deserves further investigation. 
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